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The stochastic Gross-Pitaevskii equation is shown to be an excellent model for quasi-one- 
dimensional Bose gas experiments, accurately reproducing the in situ density profiles recently ob- 
tained in the experiments of Trebbia et al. [Phys. Rev. Lett. 97, 250403 (2006)] and van Amerongen 
et al. [Phys. Rev. Lett. 100, 090402 (2008)], and the density fluctuation data reported by Armijo et 
al. [Phys. Rev. Lett. 105, 230402 (2010)]. To facilitate such agreement, we propose and implement 
a quasi-one-dimensional stochastic equation for the low-energy, axial modes, while atoms in excited 
transverse modes are treated as independent ideal Bose gases. 
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I. INTRODUCTION 

Ultracold atomic gases are proving to be extremely 
useful tools for synthesizing low dimensional quantum 
models owing to the huge degree of controlability they 
offer The effective system dimensionality may be 
tuned in experiments by manipulation of external trap- 
ping potentials, with the underlying physics ultimately 
set by the trap geometry and level of quantum degener- 
acy [2| . Increasing the trapping potential in one direction 
leads to an effectively two-dimensional (2D) system 
which allows access to a number of interesting phenom- 
ena such as the Berczinskii-Kosterlitz-Thoulcss transition 
and related studies on the nature of the Bose gas in two- 
dimensions [6T-ll3j. Increasing the trapping potential in 
a further direction results instead in an effectively one- 
dimensional (ID) system fl4l - l42j |. 

In a ID set-up, one may obtain (43j either a weakly- 
interacting system, or, for rather low densities, a 

strongly-interacting Tonks-Girardeau gas (lil - fill |44j. 
The finite temperature phase diagram of a weakly inter- 
acting ID Bose gas [H|,|43|,|4i| is more complex than that 
of a 3D gas, due to a separation in the temperatures for 
the onset of density and phase fluctuations. Density fluc- 
tuations are typically suppressed at higher temperatures 
than phase fluctuations, allowing for the formation of a 
so-called quasi-condensate [461 ] . A number of experiments 
have recently probed the physics of highly-elongated fi- 
nite temperature Bose gases at equilibrium, including the 
direct observation and analysis of density [lTl [20l. [22| and 
phase fluctuations (23T-i27j . Understanding the properties 
of matter waves in such geom etries is of key importance 
to atom interferometers j28l430l . 13214361 ] and atom chips 

In this paper, we show that in situ density profiles and 
density fluctuations from the elongated Bose gas experi- 
ments of Trebbia et al. [l8| , van Amerongen et al. [l9j , and 
Armijo et qZ.[20|. can be predicted ab initio by means of 
an effective ID stochastic model; to achieve this, we pro- 



pose and implement for the first time a modification to 
the usual form of the stochastic Gross-Pitaevskii equa- 
tion j47l - l49| which additionally incorporates quasi-lD ef- 
fects. Such an extension beyond the purely ID limit is 
required, since these experiments probe regimes for which 
both fi, ksT > huj±, as evident from FigfT] Thus, one 
should in general account for both a quasi-lD degenerate 
system exhibiting fluctuations and for the non-negligible 
role of transverse thermal modes. These are included 
here by (i) implementing a stochastic quasi-lD equation 
of state for the axial modes, and (ii) treating atoms in ex- 
cited transverse modes as independent, ideal Bose gases, 
as proposed in a related study [l9j . 

In Sec. II we discuss our methodology in more detail, 
while Sec. Ill focuses on the ab initio reproduction of the 
experimental den sity profiles and density fluctuation re- 
sults reported in jl8l - l2fjl ], with our conclusions presented 
in Sec. IV. 
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Figure 1: (Color online) Phase diagram indicating the theo- 
retically obtained parameters for the experiments considered. 
Hollow symbols indicate density profile data of [lj| and [l9l |. 
whereas filled symbols indicate the density fluctuation data 
of M- 
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II. METHODOLOGY 

We seek a model which can provide ab initio predic- 
tions for experimentally measurable properties obtained 
by in situ absorption imaging. There are two issues that 
need to be addressed simultaneously here, regarding the 
spatial extent of the quasi-condensate in the transverse 
direction, and the role of atoms in excited transverse 
modes. As mentioned, the important parameters affect- 
ing these are the ratio of the chemical potential, [i, and 
thermal energy, fcsT, to the transverse ground state en- 
ergy huij_. For fi <C huij_, the transverse ground state 
density has a Gaussian profile of width the trans- 
verse oscillator length, whereas for larger fi, this width 
becomes increased due to the effect of repulsive interac- 
tions. On the other hand, if ksT <C hu)±, then thermal 
occupation of the transverse excited modes is negligibly 
small; for higher temperatures, this is no longer true and 
atoms in these modes will contribute considerably to ex- 
perimental observables, such as density profiles. 

In the present work, we treat the quasi-condensate 
and thermal modes separately: consideration of trans- 
verse thermal modes is crucial for matching total atom 
numbers and density profiles, however they may be sim- 
ply and accurately described as independent equilibrium 
Bose gases, as shown in [l9[; fluctuating axial modes are 
instead treated within our modified stochastic model, a 
novel feature of this work, which is discussed in detail 
below. 

Quasi-condensate Modeling - The axial modes of a 
sufficiently elongated Bose gas are typically subject to 
phase and density fluctuations due to thermal excita- 
tions with wavelengths greater than the radial extent of 
the gas [13, [H| • This requires the inclusion of such fluc- 
tuations when describing modes with energies less than 
huj_. We therefore choose to describe these axial modes 
dynamically using a stochastic Gross-Pitaevskii equation 
(SGPE) [47i-[49| . In solving this equation, we make the 
following two assumptions: (i) thermal modes (with ener- 
gies > huj±) may be treated as though at equilibrium, and 
therefore represent a heat bath in contact with the axial 
sub-system (these two components are assumed to be in 
diffusive and thermal equilibrium with a temperature T 
and chemical potential /Lt); (ii) the modes in the weakly 
trapped, axial direction are sufficiently highly occu pied 
that the classical field approximation is valid [48l . |52h56| . 

Under assumptions (i) and (ii), the axial modes may 
be represented by the ID SGPE, 



inM = (l-rr(M)) 
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where tp is a complex order parameter, V(z) = mui 2 z 2 /2 
is the axial trapping potential, g = 2hu±a is the one- 
dimensional interaction strength (with a the s-wave scat- 
tering length), and rj is a complex Gaussian noise term, 



with correlations given by the relation {if{z, t)i](z' , t')} = 
2/17(2, t)k B TS(z - z')S(t - t'). The strength of the noise, 
and damping, due to contact with the transverse ther- 
mal modes, is given by 7(2,/). This may be calculated 
ab initio in terms of the Keldysh self-energy [U IH, [Hf , 
however as we are interested only in the system prop- 
erties at equilibrium, we may approximate this quantity 
to be spatially and temporally constant 85|. To a good 
approximation, this is given by 7 = 3 x 4ma 2 kBT/ (irh 2 ) 

[EE Hi- 

Eq.JT]) is valid in the scenario that the transverse 
ground state is a Gaussian of width l±. However, as we 
wish to consider experimental data from actual quasi- 
1D systems, we should additionally modify our model in 
a manner which accounts for the transverse swelling of 
the Bose gas due to repulsive interactions, as seen ex- 
perimentally [5^ . In the context of the ordinary Gross- 
Pitaevskii equation, one may replace the ID equation of 
state n[n] = gn, where n denotes the density, with [60j 



[n] = tlijj± [\/T+ Ana - l] , 



(2) 



which is obtained variationally by minimising the 3D 
GPE with respect to the transverse chemical potential 
[6lj - l63| . This was shown in l6ll t o interpolate smoothly 
across the lD-to-3D crossover [64[ , and reduces to the ID 
result in the limit that 4an <C 1, as pointed out (for the 
ordinary GPE) in [10, IH, HI]. 

Motivated by this we propose, somewhat heuristi- 
cally, a similar amendment to the ID stochastic equation 
(Eq.([T])); for elongated but not truly ID atomic clouds, 
this gives rise to the modified stochastic equation 
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which we henceforth refer to as the quasi-lD SGPE. The 
proposed modification to the ID SGPE is likely to be- 
come important when the inequality fi <C huj± is no 
longer satisfied, signaling the onset of quasi- ID effects. 
This method is of course only intended for modeling very 
elongated systems with [i < few hu±, in the weakly- 
interacting regime mg/h 2 n <C 1; in this work we thus 
restrict the application of this equation to the weakly 
interacting regime, where numerous experiments exist. 

In our stochastic scheme, the equilibrium state is 
reached in a dynamical manner, when the effects of the 
noise and damping terms of Eq.((3]) balance out. Al- 
though we demonstrate that this leads to accurate equi- 
librium predictions, the validity of this equation for de- 
scribing dynamical features remains to be investigated. 

Thermal Transverse Modes: Atoms in the transverse 
modes are considered to be in static equilibrium, dis- 
tributed according to Bose-Einstein statistics; the rela- 
tive importance of their contributions depends on the 
ratio fc^T '/hui±, and their contribution is significant in 
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most experimentally relevant cases [18M20j | (except at ex- 
tremely low temperatures). To account for their contri- 
bution to total linear density profiles, we compute (mak- 
ing the [i and T dependence explicit for clarity) 



n(z;n,T) = (\tp(z; ^,T)\ 2 ) + n±(z; (t,T) 



(4) 



where 



1 



OO 



,(M-V(z)-jHu±)/k B T 



ni _{z-^T) = —^{ 3 + l)g 1/2 

3=1 

(5) 

gi/2[- • •] is the polylogarithm (or Bose function) of order 
1/2, and AdB = hj \[2wmkBT is the thermal de Broglic 
wavelength. 

Such an addition of particles in thermal modes has 
already been used in the so-called modified Yang- Yang 
model of [l9| , where Eq. ([5]) was used in conjunction with 
a description of atoms within axial modes based on Yang- 
Yang theory. In our approach, we effectively replace the 
Yang- Yang theory with the SGPE, thereby also providing 
an indirect comparison between SGPE and Yang- Yang, 
in the weakly- interacting limit. Furthermore, on making 
this replacement, the contribution given by n± remains 
consistent with our treatment of bosons in excited trans- 
verse modes as having already reached thermal equilib- 
rium, as assumed within the SGPE model we apply here. 

To summarize briefly, our approach for modeling equi- 
librium properties of finite temperature quasi-lD Bose 
gases is based on self-consistcntly solving Eq.Q and 
Eq.([5]) for the desired total atom number (or measured 
peak density) and temperature [86j . 

In Sections IIII Al and lHI B[ we give a quantitative com- 
parison between the proposed method and published in 
situ experimental data, thereby highlighting the useful- 
ness of this method. For a direct comparison to these 
experiments, we choose here to fix the temperature to 
that reported in the experimental papers, and then use 
\i as a free parameter, which we vary until the required 
linear density is obtained at the trap centre. 



III. COMPARISON TO EXPERIMENTS 

A. Density Profiles 

We begin with a comparison of the model to total linear 
density profiles, as obtained by in situ absorption imaging 
within two experiments: Sec. IIII Al gives a comparison to 
the data of Trebbia et aZ.[l^], whose published analysis 
was based on a 3D Hartee-Fock model, before discussing 
the measurements of van Amerongen et al.fl^. who in- 
stead analyzed their results via the modified Yang- Yang 
theory. 

In each study, experimental data was compared to the- 
ory in order to quantify the equilibrium state through 



a chemical potential and temperature. The tempera- 
ture may be quite straight-forwardly measured by fit- 
ting the wings of the density distribution to the ideal 
gas result (for sufficiently high temperatures), implying 
interactions have little effect on measurements of this pa- 
rameter. Conversely, the chemical potential is far more 
dependent on the model used in analyzing the density 
profile; this is because different theories represent the full 
quantum Hamiltonian of the interacting system to differ- 
ent levels of approximation (see e.g. [65[), so incorporate 
many-body effects to a differing degree. 

It is worth pointing out also, that the density pro- 
files of the theory presented here, and those from ex- 
perimental absorption imaging share a common feature: 
namely, that additional analysis is required to identify a 
phase coherent (or 'true' condensate) and density coher- 
ent (or quasi-condensatc) fraction from the total density. 
In the SGPE, the total density is due to both coher- 
ent and incoherent particles, however knowledge of the 
first and second order correlation functions was shown to 
be sufficient to isolate both quasi-condensate and true- 
condensate densities (6(| [6?]]. In particular, the method 
discussed in (6(| may be easily applied to directly ex- 
tract such components from experimental measurements 
of correlation functions, thereby offering a more accurate, 
experimentally self-consistent, characterisation of phase- 
fluctuating experiments (without the need to resort to 
bimodal fits which become somewhat inaccurate in this 
limit). 



1. Comparison to work of Trebbia et al. \l 



The results of Trebbia et al. demonstrated experimen- 
tally the breakdown of the Hartree-Fock method when 
applied to highly elongated Bose gases [l8[. It was con- 
cluded that this breakdown occurs because density fluc- 
tuations are not accounted for accurately within this the- 
ory, since the energy lowering effect of spatial correlations 
between atoms is not captured. These density correla- 
tions [68l - l7lj are key to correctly predicting the onset of 
quasi-condensation and the associated reduction in den- 
sity fluctuations; it was found, therefore, that the excited 
states did not saturate within Hartree-Fock theory and 
so no quasi-condensate was predicted to form. 

If we compare instead to the results of the quasi- ID 
SGPE, with the total linear density given by Eq. Q, 
we see from FigJ5] that the theoretical results (black 
solid line) match well those obtained within the exper- 
iments (red circles). The agreement is extremely good 
across the entire range of temperatures considered, no- 
tably even at a temperature close to the crossover from 
quasi-condensate to thermal gas (Fig. [2£c)). It is pre- 
cisely this regime of critical fluctuations in which a de- 
viation from mean field theory might be expected, due 
to the lack of a well defined mean field quantity. Inter- 
estingly, behaviour suggestive of this was found in [l8j 
when comparing their data to the Hartree-Fock mean 
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Figure 2: (Color online) Total linear density profiles from the quasi-lD SGPE model (Eq.Q; black solid line) versus data from 
the experiment of Trebbia et aZ.[l8|] (red circles). Insets: quasi- ID SGPE contribution (Eq. (J3j ; green dashed shaded region) to 
the total linear density profiles (Eq.([4|; black solid line). Bottom: Table showing parameters /i and T for the quasi-lD SGPE 
density profiles which match the experimental data; the parameters given in from a 3D Hartree-Fock fit are also shown. 



field model: the experimental data was found to have 
a higher peak than the mean field result [see Figl(c) of 
[18|]. which illustrates the potential importance of includ- 
ing many-body effects, as studied recently in the context 
of a finite-temperature classical field theory 67 1. 

We obtain a good fit between the quasi-lD SGPE 
and experimental density profiles, at the experimentally 
measured temperatures, with a comparison between the 
chemical potentials extracted in our treatment and the 
published values based on the Hartree-Fock analysis of 
[18j shown in the table |87[- The deviation between the 
parameters of different theoretical methods should not 
be of any concern, as it merely highlights fi as a model 
dependent quantity, i.e. it is dependent on the actual 
Hamiltonian used to analyse the experimental results, 
and therefore varies depending upon the level of approx- 
imation 1661 16711 . 



The insets of Fig(5] show the contributions to the to- 
tal linear density profiles due to the axial SGPE density 
(green dashed, shaded region), with the remainder com- 
ing from Eq.©. The importance of the SGPE contribu- 
tion is clear in the first three plots, which shows an ap- 
preciable fraction of atoms reside in axial modes for these 



parameters. In the highest temperature case, shown in 
FigHJd), /i < and the gas is entirely in the thermal 
phase; here the density is instead due almost entirely to 
atoms in tranverse modes. 



2. Comparison to work of van Amerongen et al. flf\l 

The second experiment that we consider is that of van 
Amerongen et al. This was the first experimen- 

tal comparison to the exact Yang- Yang thermodynamic 
solution to the finite temperature ID Bose gas problem 
[72j, also referred to as the thermodynamic Bethe ansatz. 

In [l9j . the one-dimensional Yang- Yang theory was 
used to represent the axial modes and transverse ground 
state (of width while the contribution to the linear 
density due to atoms in transverse excited states, was 
accounted for using the method we also adopt here. The 
total density profiles in were therefore calculated us- 
ing Eq.©, with the role of the SGPE contribution (l^l 2 ) 
instead played by the ID Yang- Yang prediction. There- 
fore, in comparing to this work, we will gain insight on 
two fronts: firstly, how well the SGPE matches the ex- 
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Experiment of van Amerongen et al. [19] 




200 400 













Analysis of [19] 


Current analysis 




modified Yang- Yang 


ID SGPE 


Quasi-ID SGPE 


Fig. 3a 


4.0 


-2.74 


-2.74 


-2.74 


Fig. 3b 


2.8 


-0.38 


-0.38 


-0.38 


Fig. 3c 


2.4 








0.03 


Fig. 3d 


0.89 


0.45 


0.45 


0.39 



Figure 3: (Color online) Top: Total linear density profiles from the quasi-lD SGPE (black solid line) versus the experimental 
data of van Amerongen et al. pj| (red circles). Insets: Quasi-ID SGPE (thick black line) versus ID SGPE (thin brown line) 
density profiles. Bottom: Table showing parameters /i and T for the ID and quasi-lD SGPE density profiles which match the 
experimental data and modified Yang- Yang model fits from [3] . 



pcrimcntal data in this regime, and simultaneously (but 
indirectly), how well the SGPE prediction for the density 
profiles matches that due to Yang- Yang thermodynamics. 

Figj3] shows that the agreement between the experi- 
mental data and the proposed quasi-lD SGPE approach 
is again very good across the entire temperature range 
probed, including the crossover from quasi-condensate to 
degenerate thermal gas. 

We now wish to discuss how our quasi- ID SGPE re- 
sults compare to those from the SGPE with the usual ID 
equation of state. Practically, this means using the equi- 
librium result of Eq. ([T]), rather than Eq. (|3|), as the axial 
density input (IV'I 2 ) in Eq. (j4]). Although this approach 
also recovers closely the total density profiles found with 
the quasi-lD SGPE at each temperature (insets to Fig. 
|3]) , and therefore also those measured experimentally, it is 
important to note that each approach can lead to slightly 
different chemical potentials for the same temperature. 

Importantly, we find that the parameters used to ob- 
tain the ID SGPE results are identical to those obtained 
from fits of the modified Yang- Yang model to the den- 
sity data in [l9|]. These results have also been reported 
[73| to arise within the context of the closely-related ID 



stochastic projected Gross-Pitaevskii equation (SPGPE) 
[7^ in parallel independent work, which also looked at 
the momentum distribution of the gas after focussing 
[75j . This provides an indirect additional test between 
the ID SGPE, the ID SPGPE and Yang- Yang theories 
in the weakly-interacting regime. The parameters pre- 
dicted by the quasi-lD SGPE, ID SGPE and modified 
Yang- Yang models are shown in the Table of FigJ3J 

Having established that both the ID and quasi-lD 
SGPE approaches accurately reproduce experimental 
density profiles (and therefore also those due to the mod- 
ified Yang- Yang approach used in [lH), we now turn to 
an investigation of density fluctuations, which provide a 
more sensitive probe for the validity of these theories. 



B. Density fluctuations 

Density fluctuations are increased markedly within an 
ideal Bose gas, due to an effect of quantum statistics, 
which leads to atomic bunching 76]. However, at suf- 
ficiently low temperatures, and in the presence of inter- 
actions, quasi-condensation leads to atom-atom correla- 
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tions which overcome the tendency for bosonic atoms to 
bunch together, and therefore to a reduction in the level 
of density fluctuations, relative to those expected in an 
ideal Bose gas [22j . 



Comparison to work of Armijo et al. \2(\l 

In a recent paper [20] , Armijo et al. measured the sec- 
ond and third moments of the density fluctuations of a 
finite temperature Bose gas, comparing these to theoreti- 
cal predictions from ideal Bose gas and quasi-condensate 
mean-field models, and also the modified Yang- Yang 
model of . We now briefly outline their experimental 
method before describing the numerical scheme we follow 
in order to closely mimic this. 

Experimental procedure: A harmonic trapping poten- 
tial leads to a density profile in which the number of 
atoms varies spatially, and therefore close to the trap 
centre it is possible to have a scenario where there is suf- 
ficient levels of degeneracy such that a quasi-condensate 
is formed, whereas in the low-density wings, the gas 
is still effectively a non-interacting thermal gas. Thus, 
at a single temperature, by scanning the spatial extent 
of the trapped gas, it is possible to observe both the 
enhancement of density fluctuations, due to quantum 
statistics (low density, ideal Bose gas), and their sub- 
sequent suppression, due to particle interactions (higher 
density, quasi-condensate regime). 

An approach based on this observation was first under- 
taken experimentally in [22j . and subsequently followed 
by more detailed studies in (2(| In (20[, the gas 

was probed using a CCD camera, which effectively di- 
vides observations of the gas into pixel sized regions (of 
size A = 4.5pm in this case). Absorption imaging al- 
lowed for the number of atoms within each pixel, N , to 
be measured, and repeated measurements provided a set 
of fluctuating values, as well as an average number per 
pixel, (N). The p-th moment of the density fluctuations 
for the set of measurements were then calculated for each 
pixel as usual via (SN?) = ((N — (N)) p ). 

The aim of the present work is to demonstrate the 
SGPE as an ah initio method for analyzing experimen- 
tal findings, and so we wish to implement a numerical 
scheme which follows experimental procedures as closely 
as possible. Fortunately, the grand canonical formulation 
of the SGPE makes it relatively simple to simulate the ex- 
perimental methods used in [20[ . This is because, in addi- 
tion to the unified treatment of both (quasi-)condensate 
and thermal atoms in density profiles, the SGPE shares 
a second feature in common with experiments, namely 
a shot-to-shot variation between individual realizations. 
This variation enables us to straightforwardly model the 
equilibrium density fluctuation experiments of [2(j, but 
is also important in dynamical studies (77l - l80| . 

Numerical SGPE procedure: Physical observables 
within the SGPE arc obtained as products of the stochas- 
tic wavefunction ip, averaged over many realizations of 



the noise r\ (see Eq.[3l and Ref. [81| for a simple overview). 
This leads naturally to a shot-to-shot variation between 
numerical realizations, in a way analogous to an indi- 
vidual experimental run. By treating each noise real- 
ization like an experimental realization, we are able to 
accurately synthesize the experimental procedure of [20| 
within our numerical simulations. We emphasise, how- 
ever, that within this analogy it should be understood 
that single runs represent weighted contributions to av- 
erages over fluctuating quantities, and that it is only 
such appropriately averaged quantities that we expect to 
match well those determined experimentally, as indeed is 
found to be the case below. 

Our numerical procedure is to run a large number of 
stochastic simulations (1000), each of which yields a fluc- 
tuating density profile. We measure the number fluctua- 
tions within each pixel sized region, by first spatially bin- 
ning the SGPE density data into A-sized regions |88[ , in 
order to give an output consistent with that obtained in 
the experiments (see FigfJ] (a)-(b)). Integrating over the 
numerical grid points within a single pixel yields iV z a s 
the value for the (fluctuating) atom number from a sin- 
gle stochastic realization, within that pixel. Repeating 
the same procedure for the mean total density generates 
a binned average axial pixel number, N z ^. The contri- 
bution from the axial modes to the second (p = 2) and 
third (p = 3) moments of the density fluctuations are 
then calculated as (5NP) Z = ((iV Z)A - N zA ) p ). 

Similarly, a binned transverse contribution to the av- 
erage pixel number -/Vi a may be obtained, and the to- 
tal average atom number in each pixel is then given 
by Na = N Z ,A + ^V_l,a (which corresponds to (N) in 
[201). As we treat the atoms in the transverse modes 
in a static way, they give a non-zero contribution only 
to average properties, and so do not contribute to mo- 
ments of the density fluctuations directly. However, as 
we found in Section IIII Al atoms in these modes are 
well approximated as a de gen erate ideal gas, for which 
(SN 2 )± ~ (SN 3 )± ~ (N)± [20I, and we therefore assume 
that these atoms contribute a factor (N)± to both sec- 
ond and third moments. So, ultimately, we compute the 
total density fluctuations as (5N P ) = (5N P ) Z + {N)x, for 
p = 2,3. 

To illustrate our procedure, we plot in Figf?] an ex- 
ample single run and average density profile both before 
(FiglUa)), and after (FigQJb)), spatial binning. Com- 
paring these plots, it is clear that the binning procedure 
significantly smooths the raw single run data, as would be 
expected for this kind of a spatial coarse graining proce- 
dure; note that the binned data looks remarkably similar 
to the experimental fluctuating density profile shown in 
[2(| (Fig. 1(c) of that work). We additionally show in 
FiglUc) the variance in the density fluctuations which 
results from the set of 1000 fluctuating binned densities; 
this is plotted against the average number of atoms per 
pixel [89l | . 

Making use of the thermodynamic relation (SN 2 ) = 
kBTA(dn/d[i)T (zlj, it is possible to derive mean field 
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Figure 4: (Color online) Top row: quasi-lD SGPE den- 
sity (noisy red curve) obtained from (a) a single numerical 
run (raw data), and (b) corresponding spatially binned data; 
dashed black curves display the density averaged over 1000 
independent realisations, (c) {SN 2 ) m [8(| from the binned 
quasi-lD SGPE data (black diamonds) and ID SGPE data 
(green squares), compared against correponding mean field 
results for an ideal Bose gas (brown solid line) ID quasi- 
condensate (dotted maroon, horizontal) and quasi-lD quasi- 
condensate (dashed red); the thin, vertical dashed lines indi- 
cate the 'crossover' region where the interaction and thermal 
energies become comparable. This data is for T = 96nK. 



results for the density fluctuations, based on both the 
ideal gas and quasi-condensate equations of state. The 
ideal gas result is 



{SN *> = _L y VLej^ 1 

^ ' XdB^k B T V7 (l-e-^/ fc BT)' 



(6) 



while using the ID equation of state for a quasi- 
condensate (fJi[n] = gn) gives the simple result [22[ 
(<5iV 2 )i D = k B TA/g. Using instead the quasi- ID 
equation of state, Eq.J2]), yields, (<5A 2 ) quasi _ 1D = 

(<W 2 ) 1D [l + (f,-V{z))/taj ± ]. 

Within the higher density region of Fig|3tc) , we show 
the mean field results due to the ID and quasi- ID equa- 
tions of state for a quasi-condensate. It is clear that the 
ID SGPE shows good agreement with the ID mean held 
result, whereas the quasi- ID SGPE instead agrees very 
well with the quasi- ID mean held prediction. The two 
vertical lines indicate the region where the interaction 
energy and the average thermal energy become compa- 
rable. For higher densities, interactions significantly re- 
duce (5N 2 ) below the ideal gas prediction, and we see 
this is more pronounced in the ID SGPE case relative to 
the quasi-lD SGPE data, as expected from the ID and 
quasi-lD mean field predictions. 

A key point from our analysis at this temperature, is 
that while the ideal gas equation of state is valid only 
for small densities, and the mean-held quasi-lD equation 
of state holds at high densities, the quasi-lD SGPE, like 
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Figure 5: (Color online) Second (top row: (a)-(b)) and third 
(bottom row: (c)-(d)) moments of the atom number fluctu- 
ations from the quasi-lD SGPE data (black diamonds) and 
experimental data from the paper of Armijo et al. [53] (red 
circles). Temperatures are T = 376nK (left images) and 
T = 96nK (right images). Insets: axial (light blue crosses) 
versus axial plus transverse (black diamonds) contributions to 
the number fluctuations vs. the ideal gas result (brown solid 
line). 



the experimental data, provides a smooth crossover be- 
tween each of these regimes. This is because, in moving 
outwards from the centre of the trap in the presence of 
a quasi-condensate, at some point the gas changes phase 
to a thermal gas, and a mean field theory cannot be ex- 
pected to accurately describe fluctuations in the transi- 
tion region near the edge of the quasi-condensate. 

Comparison to Experiment: Fig JS] shows a comparison 
of the experimental density fluctuation data obtained by 
Armijo et al. (red circles), and the quasi-lD SGPE model 
(black diamonds). We plot in the top row (SN 2 ) m , and 
in the bottom row (6N 3 ) m , each versus for two tem- 
peratures (T = 376nK (left) and 96nK (right)). 

Concentrating on the (5N 2 ) m data first, it is clear that 
the quasi-lD SGPE numerical results follow the exper- 
imental data well in both the high temperature (left) 
and low temperature cases (right images). Notice that 
as the high-density, quasi-condensate regime is reached 
in FigEJb), the ID SGPE result (shown in FiggJ} would 
predict too great a reduction in density fluctuations rela- 
tive to the experimental results, whereas the quasi-lD 
SGPE captures the experimental behaviour very well. 
Physically, this suggests that the effect of the transverse 
swelling of the quasi-condensate near the centre of the 
trap cannot be ignored for these parameters, and that 
the quasi-lD extension to the SGPE is therefore essen- 
tial here. 

Armijo et al. found a similar trend when comparing 
between their experimental data and the modified Yang- 
Yang model. While the Yang- Yang result matched the 
low density data well, at higher densities it displayed the 
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same behaviour as the ID SGPE result. Instead, Armijo 
et at found the quasi-lD mean field result to better 
reproduce the experimental behaviour in the higher den- 
sity quasi-condensate regime. So, as for the density pro- 
files, we again find the ID SGPE data to reproduce the 
ID Yang- Yang behaviour, whereas the quasi-lD SGPE 
goes beyond the modified Yang- Yang model by captur- 
ing the quasi-lD behaviour found in the experiment in 
all regimes probed. 

The results for the third moment (FigO bottom row), 
(SN 3 ) m , also show good agreement between the quasi- 
lD SGPE and the experimental data (within the large 
experimental error bars). 

The insets to Figs. [S][a)-(b) show clearly that the re- 
duction in density fluctuations compared to the ideal 
gas prediction (brown solid line) is more pronounced in 
the low temperature case, where there is a higher quasi- 
condensate fraction. Contrary to this, the enhanced pres- 
ence of thermal atoms in the high temperature case, 
provides an opportunity to test whether our approxi- 
mate treatment of the transverse mode contribution to 
(SN 2 ) m , leads to the correct behaviour at low density. 
The results before (light blue crosses) and after (black di- 
amonds) the transverse contribution is added are shown 
in the insets to FigsG^a)-(b). This addition results in 
a noticeable shift upwards in the data of FigJSJa) while 
for the lower temperature case, the difference following 
this addition is negligible. Comparing the SGPE data 
to the expected low density result, the ideal gas predic- 
tion for the same parameters (solid brown line), we see 
that adding a contribution (N) j_ to the axial contribution 
fully captures the expected behaviour. 

IV. CONCLUSIONS 

In conclusion, we have proposed and implemented 
a suitably modified one-dimensional stochastic Gross- 
Pitaevskii equation, which was shown to provide excel- 
lent ab initio predictions for both in situ experimen- 
tal density profiles obtained by Trebbia et aZ.[l8[ and 
van Amerongen et al.fisty. and in situ density fluctuation 
data from the experiment of Armijo et aZ.[20|. This was 



achieved by matching peak densities (equivalent to to- 
tal atom number) to a hybrid scheme, which combines 
the aforementioned stochastic model with a previously 
reported approach based on treating transverse thermal 
modes as independent ideal Bose gases. 

The study of density fluctuations showed that our com- 
bined approach captures all experimental regimes stud- 
ied in a unified manner, smoothly interpolating between 
mean field models, whose individual validity is restricted 
to either the low density or high density regimes. Impor- 
tantly, it was found that analyzing individual stochastic 
realizations in the same way as individual experimental 
runs, led to good agreement between the density statis- 
tics in each case. 

Reducing our stochastic model to the previously 
tested one-dimensional stochastic Gross-Pitaevskii equa- 
tion showed that: (i) the latter model is consistent with 
Yang- Yang predictions (in the weakly-interacting regime 
probed here), and that (ii) while both one-dimensional 
and quasi-one-dimensional approaches accurately repro- 
duce equilibrium density profiles, they do so with slightly 
different chemical potentials. 

The confidence gained from these successful compar- 
isons suggests that the quasi-one-dimensional stochastic 
Gross-Pitaevskii model, which was proposed here as a 
hybrid of different previously implemented approaches, 
is an excellent model for describing quasi-condensate ex- 
periments at equilibrium; the extent to which this model 
can be directly applied to study dynamical features, such 
as experimentally-measured properties after expansion, 
remains to be investigated. 
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